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The static Casimir effect describes an attractive force between two conducting plates, due to 
quantum fluctuations of the electromagnetic (EM) field in the intervening space. Thermal fluctu- 
ations of correlated fluids (such as critical mixtures, super-fluids, liquid crystals, or electrolytes) 
are also modified by the boundaries, resulting in finite-size corrections at criticality, and additional 
forces that effect wetting and layering phenomena. Modified fluctuations of the EM field can also 
account for the 'van der Waals' interaction between conducting spheres, and have analogs in the 
fluctuation-induced interactions between inclusions on a membrane. We employ a path integral 
formalism to study these phenomena for boundaries of arbitrary shape. This allows us to examine 
the many unexpected phenomena of the dynamic Casimir effect due to moving boundaries. With 
the inclusion of quantum fluctuations, the EM vacuum behaves essentially as a complex fluid, and 
modifies the motion of objects through it. In particular, from the mechanical response function 
of the EM vacuum, we extract a plethora of interesting results, the most notable being: (i) The 
effective mass of a plate depends on its shape, and becomes anisotropic, (ii) There is dissipation 
and damping of the motion, again dependent upon shape and direction of motion, due to emission 
of photons, (iii) There is a continuous spectrum of resonant cavity modes that can be excited by 
the motion of the (neutral) boundaries. 



I. OUTLINE 



Fluctuation-induced forces are ubiquitous in nature, 
covering many topics from biophysics to cosmology . 
There are two basic ingredients in these phenomena: 
(i) A fluctuating medium, such as the electromagnetic 
(EM) field; and (ii) External objects whose presence sup- 
presses (or in some way modifies) the fluctuations, such 
as dipoles or conductors. The overall strength of the in- 
teraction is proportional to the driving energy of fluctua- 
tions (/csT and h for thermal and quantum fluctuations, 
respectively); its range is related to that of the correla- 
tions of the fluctuations. The most interesting cases are 
when the interactions are long-ranged, corresponding to 
scale free fluctuations. 

The goal of this article is to provide a glimpse of the 
unity and simplicity of fluctuation-induced forces. While 
we attempt to describe a wide range of phenomena, this 
selection is by no means exhaustive, and highly biased by 
subjective interests. In the spirit of a colloquium, we have 
tried to avoid technical details, preferring to present di- 
mensional arguments whenever possible. The interested 
reader is referred to various sources for calculational de- 
tails. 

A prototype of fluctuation-induced interactions is the 
Casimir force between conducting plates, due to quan- 
tum fluctuations of the EM fleld. In Secjl| we discuss 
several cases were the source of interaction is the thermal 
fluctuations of a correlated fluid between the bounding 
plates. This interaction was originally proposed for a liq- 
uid mixture at its critical point, but is also present when 



long-range correlations appear as a result of symmetry 
breaking, as in superfluids or liquid crystals. There is 
even an attractive component to the (mainly repulsive) 
force between two similarly charged plates, due to fluctu- 
ations of counterions in a neutralizing solution. Since the 
latter connection is seldom made explicit, we expand on 
its origin in Appendix While the reader may skip any- 
one of these sections without losing general track of the 
article, we note that experiments on wetting of helium 
fllms may provide a beautiful test of these forces. 

The van der Waals and London dispersion forces be- 
tween atoms and molecules can also be attributed to the 
modifled fluctuations of the EM fleld. As we point out in 
Sec. Ill, there are analogous forces between inclusions on 
a membrane. Their origin is the modifled surface fluctu- 
ations, and they decay more slowly with separation than 
the standard van der Waals interaction. We use this ex- 
ample to emphasize that non-additivity is an important 
feature of fluctuation-induced forces: they cannot be ob- 
tained from a pairwise sum of two-body potentials. 

Many new results are obtained in going beyond the 
simple geometries of flat plates and simple spheres, 
by looking at rough and deformed structures, as in 
Sec. IV. Our key to implementing the corresponding 



non-standard boundary conditions is a path integral ap- 
proach, which is sketched in Appendix This approach 
can also be used in conjunction with a path integral quan- 
tization of the electromagnetic field. Since in this rela- 
tivistic theory, deformations in space and time appear on 
the same footing, we can examine the dynamic Casimir 
effect which is introduced in Sec.^. 
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Some of the unexpected phenomena that emerge from 
quantum fluctuations of the EM fleld in the presence of 
moving deformed plates are discussed in Sec.VI. There 
are corrections to the mass of a plate that depend on 
its shape. There is also dissipation due to emission of 
photons (which accounts for the 'friction' in the title of 
this article). While these effects are typically very small, 
we believe that they are significant for what they imply 
about the nature of the quantized EM vacuum. Quali- 
tatively, with the inclusion of quantum fluctuation, the 
vacuum behaves as a complex fluid which influences the 
bodies moving through it. 



II. FLUCTUATION-INDUCED FORCES 

A. Quantum fluctuations 

The standard Casimir effect is a macroscopic man- 
ifestation of quantum fluctuations of vacuum. In 1948, 
Casimir considered the electromagnetic field in the cavity 
formed by two conducting plates at a separation H . Be- 
cause the electric field must vanish at the boundaries, the 
normal modes of the cavity are characterized by wave- 
vectors k = {kxjkyjTm/H), with integer n. Once quan- 
tized, these normal modes can be regarded as harmonic 
oscillators of frequencies u^k) = c\k |; each of which in its 
ground state has energy huj{k)/2. While adding up all 
the ground state energies leads to an infinite contribution 
to the overall energy £{H), Casimir showed that a finite 
attractive force is obtained from 



F(H) = -- — = -he X — 7 X , 

^ ' dH 240' 



(1) 



where A is the area of the plates. Thus, by measuring 
the mechanical force between macroscopic bodies, it is in 
principle possible to gain information about the behavior 
of the quantum vacuum. 

The predictions of Casimir where followed by experi- 
ments on quartz and aluminum |^ plates at separa- 
tions H > lO'^A. However, these experiments, and others 
reviewed in Ref. provided results that were at best 
in qualitative agreement with Eq.(|l|). It is only quite re- 
cently that high precision measurements of the force (us- 
ing a torsion pendulum) between a gold plate, and gold 
plated sphere, confirmed the accuracy of the theoretical 
prediction to within %5 M. 



B. Thermal fluctuations 

While the Casimir interaction is due to the quantum 
fluctuations of the electromagnetic field, there are several 
examples in classical statistical mechanics, where forces 
are induced by the thermal fluctuations of a correlated 
fluid. One of the best known examples comes from the 



flnite size corrections to the free energy at a critical point 
[Q . Fisher and de Gennes [|lO| argued that in a binary liq- 
uid mixture, the concentrations near a wall are perturbed 
only over a distance of the order of the correlation length 
^. Any interaction mediated by the concentration fluc- 
tuations must also decay with this characteristic length. 
However, at the critical point, where ^ diverges, they 
suggested an attractive contribution to the free energy of 
a critical fllm, that varies with its thickness iJ, as 



ST{H) = -ksT X A X A. 



(2) 



This is to be expected on dimensional grounds, as the 
free energy comes from thermal fluctuations, hence pro- 
portional to /csT, and must be extensive in A. (Sim- 
ilar analysis in d-dimensions leads to a dependence as 
1/H'^~-^.) In two dimensions, exact values for the di- 
mensionless amplitude A can be obtained by employing 
techniques of conformal fleld theories ||Tl[| . In higher di- 
mensions, they can be estimated numerically [ p2[ , and by 
e = A ~ d expansions . 

In analogy to the Casimir energy, we can regard Eq. (|^) 
as due to the modifled free energy of concentration fluctu- 
ations by the boundaries. However, the force that results 
from this free energy decays as The difference 

in power of H from Eq.(|l|) is explained by noting that 
the fluctuation energy in the latter is quantum in origin, 
hence proportional to he, which has dimensions of energy 
times length. 



C. Superfluid films 

In fact, long-range forces are induced by thermal fluc- 
tuations of any correlated medium, by which we mean 
any system with fluctuations that have long-range corre- 
lations The critical system is a very particular ex- 
ample; much more common are cases where long-range 
correlations exist due to Goldstone modes of a broken 
continuous symmetry, as in superfluids or liquid crystals. 
A superfluid is characterized by a complex order param- 
eter, whose phase 4> may vary across the system. The 
energy cost of such variations is governed by the Hamil- 
tonian 



K 



j d3x(V0)' 



(3) 



where the "phonon stiffness" K is related to the super- 
fluid density. There is no characteristic length scale for 
fluctuations of (/>, which scale as a power of the obser- 
vation length. Consequently, we expect power law finite 
size scaling, just as in the case of a critical point. In the 
Casimir geometry, the free energy resulting from thermal 
fiuctuations of these modes has the form iistl 



(4) 
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Note that the result is universal, i.e. independent of 
the stiffness K. A similar expression is obtained for the 
free energy of the electromagnetic field confined between 
metallic plates at high temperatures fcsT ^ Tic/H. How- 
ever, the result is larger by a factor of two |]l6|, reflecting 
the two polarizations of the normal modes (photons). 

Liquid Helium is a powerful wetting agent, which tends 
to spread over most surfaces. The thickness of the wet- 
ting layer is controlled by the strength of the attractive 
forces that bind the fllm to the substrate [0 , mostly due 
to van der Waals interactions. In the presence of a chem- 
ical potential penalty of Sfj, per unit volume, the energy 
of a film of thickness H is 



E{H) 



= A 



H 



(5) 



where C is a positive numerical constant. Minimizing 
this expression leads to a thickness jlj] 



„ (2CkBT 



1/3 



(6) 



When the helium film is in the normal phase, the 
film thickness is determined solely by the strength of 
the van der Waals (vdW) force. The numerical value 
of C> = Cvdw > depends on the substrate, and 
is nonuniversal. However, when the film becomes su- 
perfluid, there is an additional attractive fluctuation- 
induced (FI) force due to Eq.(|4|), and C< — Cydw + Cfi', 
where Cpi = -C(3)/167r « -0.02391. In the vicin- 
ity of the superfluid transition, there is a different at- 
tractive contribution to the force due to finite size scal- 
ing (FSS) of the critical fluctuations, as in Eq.(^), and 
C\ — Cvdw + Cfss- The best estimate for the finite 
size scaling amplitude for the XY model in d = 3 is 
Cfss ~ —0.03 |l^,^. The parameter C thus takes three 
different values in the normal fluid, at the A-point, and 
in the superfluid phase. From Eq.(|^) we then expect two 
jumps in the film thickness, as the temperature is low- 
ered through the superfluid transition. Experiments to 
monitor the film thickness, thus providing a test of these 
forces, are currently underway at Penn. State University 



D. Liquid crystals 

Liquid crystals exemplify anisotropic cases of corre- 
lated fluids due to broken symmetry, which again lead to 
fluctuation-induced forces ||2^,^,|l5|. They are also eas- 
ily accessible, as experiments can be performed at room 
temperature and require no fine tuning to achieve crit- 
icality. A nematic liquid crystal is composed of long 
molecules that are aligned, with an order parameter 
which is the 'director' field n(r), characterizing the lo- 
cal preferred direction of the long axis of the molecules 
p2| . The energy cost of fluctuations of this field is given 
by [H 



N = ^Jd^r[ Ki(V • n)2 + K2(n • V x n)^ 

-KK3(n X V X n)2]. (7) 



Integrating over the nematic fluctuations leads to a free 
energy contribution 



5£n = -ksT X ^ X ^ 



Hi 



11 

K2 



(8) 



Note that the resulting force does depend on the relative 
strengths of the elastic coupling constants (reflecting the 
anisotropy of the system) . 

In a smectic liquid crystal, the molecules segregate into 
layers which are fluid like. The fluctuations of these lay- 
ers from perfect stacking are described by a scalar defor- 
mation m(x, z), which is subject to a Hamiltonian 



oz 



The resulting interaction energy 

S£s = -kBT X A X 

MX Idtt 



with A 



(9) 



(10) 



falls off as 1/H, reflecting the extreme anisotropy which 
has introduced an additional length scale A into the prob- 
lem. For potential experimental tests of these forces in 
surface freezing of liquid crystal films, see Ref. [ESl. 



E. Charged fluids 

Interactions between a collection of charged macroions 
in an aqueous solution of neutralizing counterions, with 
or without added salt, are in general very complex. 
The macroions may be charged spherical colloidal par- 
ticles, charged amphiphilic membranes, stiff polyelec- 
trolytes (e.g. microtubules, actin filaments, and DNA), 
or flexible polyelectrolytes (e.g. polystyrene sulphonate), 
and the counterions could be mono- or polyvalent. It 
is known that, under certain conditions, the accumula- 
tion (condensation) of counterions around highly charged 
macroions can turn the repulsive Coulomb interaction be- 
tween them into an attractive one. The attractive inter- 
action is induced by the diminished charge-fluctuations 
close to the macroions, (due to the condensation of coun- 
terions) |p^-p8||, and in this sense related to the effects 
discussed in the previous sections. 

Since the connection between the entropic attraction of 
charged macroions and the general class of fluctuation- 
induced forces, is seldom made explicit, in Appendix 
we present a path integral formulation that makes 
this analogy more transparent. The interaction between 
macroions can be broken into two parts: A Poisson- 
Boltzmann (PB) free energy, and a fluctuation-induced 
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correction. Specifically, consider two parallel negatively 
charged 2D plates with densities — cr, separated by a dis- 
tance 77 in d = 3, in a solution of neutralizing counterions 
with valence z. The PB equation can be solved exactly 
in this geometry, and the corresponding PB free energy, 
in the limit of highly charged plates, is 12^ 



A 



PB 



1 



1 



(11) 



in which is = /eksT is the Bjerrum length. Note that 
in the limit iszaH 3> 1, the interaction is independent 
of the charge densities of the plates; i.e. it is universal. 

The fluctuation-induced correction involves calculation 
of a determinant (see Appendix ||), which depends on 
the local charge compressibilities. The true compress- 
ibility profile (and the charge density profile) emerging 
from of the solution of the PB equation, is generally 
very complicated. It is usual to simplify the problem 
by assuming that the surface charge density is so high 
that the counterions are confined to a layer of thickness 
A <C -ff , where A ^ 1/z£b<^ is the Gouy-Chapman length. 
Then we can use an approximate compressibility pro- 
file m2(z) = {2/0[S{z + H/2) + 5{z - H/2)], in which 
= TT^i'^z^cr^A defines a "crossover length". In the 
limit H/^ 1, we again obtain 



Ffi = -ksT X A X ^ [1 + 0{^/H)] , (12) 



for the fluctuation-induced part of the interaction plj 



III. DISPERSION FORCES 
A. Van der Waals interactions 

In addition to his work on the force between plates, 
Casimir also realized that the van der Waals and 
London |33j forces can also be understood on the same 
footing: The presence of the atoms modifles the fluctu- 
ations of the electromagnetic field, resulting in an at- 
tractive interaction. (For a modern perspective, see the 
discussion by Kleppner in Ref. |Q.) 

For example, let us consider two conducting spheres of 
volumes Vi = 47raf/3 and V2 = 47ra2/3, at a distance 
R. Usually, the van der Waals interaction is obtained 
from the instantaneous induced dipoles on the spheres. 
However, it can be equivalently obtained by examining 
the electromagnetic fluctuations of the remaining space. 
The fluctuation induced interaction is proportional to the 
product of the excluded volumes, and thus on dimen- 
sional grounds we expect a potential 



V{R) = -keT X X At, 



(13) 



due to thermal fluctuations. When the fluctuations are 
of quantum origin, Eq.(|l^) is modiflcd to 



V^V■? 

V{R) = -hcx^xA^ 



(14) 



with Aq = 23/(47r) H. 

Let us compare the above result with the more stan- 
dard approach to calculating the London force between 
two neutral atoms: While the average dipole for each 
atom is zero, an instantaneous dipole fluctuation in one 
can induce a parallel instantaneous dipole on the other, 
leading to an attraction. Since the direct dipole-dipole 
interaction decays as 1/i?^, the induced effect scales as 
the square, i.e. 1/R^. In this regards, it is similar to the 
result in Eq. ([l^) , except that the characteristic energy is 
set by a typical atomic excitation energy of huo rather 
than ksT. The retardation effects are then obtained by 
taking into account the finite speed of light. We can 
imagine that for large enough distances, when the sig- 
nal goes from atom-1 to atom-2 to induce the dipole, 
and return to atom-1 to induce an attraction, it finds 
the dipole at atom-1 somewhat misaligned; resulting in 
a weaker attraction. The characteristic time for electron 
movements can be estimated from the frequency of the 
orbit as T = I-k jio^. The crossover occurs when the travel 
time for the signal is comparable to this characteristic 
time, namely Rjc ^ t. Hence, taking into account the 
retardation effect, the interaction is 



Vr{R) = -hivo X X / ( — 



(15) 



The crossover function f(x) is a constant for x — > 0, and 
vanishes as 1/x for a; 00. In the latter limit, the de- 
pendence on LOo vanishes, and the Casimir-Polder result 
of Eq. ( p^ ) is recovered. A similar crossover function be- 
tween the two forms of interaction for conducting spheres 
in Eqs.(p^) and (p^, occurs at a distance At ~ hc/kBT. 



B. Inclusions on membranes 



Dispersion forces are not limited to particles in three 
dimensional space, but also occur for inclusions on films 
and membranes, the latter is of potential importance for 
understanding the interactions between proteins fioating 
on a cell membrane. A membrane is a bilayer of am- 
phiphilic molecules, each composed of a hydrophilic or 
polar head, and a hydrophobic tail of hydrocarbon chains. 
The polar heads prefer to be in contact with the water, 
and in the bilayer structure insulate the "oily" hydrocar- 
bon interior from contact with water. A bilayer that is 
in equilibrium with amphiphiles in solution, can easily 
change its area by exchanging molecules with this reser- 
voir. This implies that the surface tension is zero [pO|-^, 
and the energy cost of deforming the bilayer is entirely 
due to bending |39|. The "fiicker" of the membrane is 
thus governed by the elastic Hamiltonian 
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(16) 



where ft.(x) is a height function describing deformations 
of the surface. Typical values of the bending rigidity k, 
of biological membranes is of the order of 10^ °K. 

Cell membranes also include proteins that perform var- 
ious biological functions (e.g. pumps). Each protein in- 
clusion disturbs the lipid bilayer, resulting in interactions 
between nearby inclusions (c.f. Refs. [p5|-p8[ and ref- 
erences therein). These disturbances, and the resulting 
interactions, tend to be short-ranged, falling off exponen- 
tially with a characteristic length related to the distance 
over which the lipid membrane "heals" |Q. There are 
also longer ranged interactions between the proteins: In 
addition to the standard van der Waals interaction, there 
are interactions mediated by the disturbed fluctuations of 
the flickering membrane. As discussed in Ref. |^^, such 
interactions exist as long as the rigidity of the inclusion 
differs from that of the ambient membrane, and fall off 
as l/R^. In particular, if the inclusions are much stiffer 



than the membrane, the fluctuation-induced potential is 



V(i?) 



6 



(17) 



where A is the area of each inclusion Q] . The interaction 
is attractive and independent of k and k; its energy scale 
set by ksT. A generalization of this result that includes 
quantum fluctuations of membranes is given in Ref. . 

The form of the interaction depends sensitively on the 
shapes of the inclusions, as demonstrated by the calcula- 
tion of the fluctuation-induced interaction between rod- 
like objects The rods are assumed to be sufhciently 
rigid so that they do not deform coherently with the un- 
derlying membrane. They can thus only perform rigid 
translations and rotations while remaining attached to 
the surface. As a result, the fluctuations of the membrane 
are constrained, having to vanish at the boundaries of the 
rods. Consider the situation depicted in Fig. |^, with two 
rods of lengths Li and L2 at a separation R ^ Li. The 
fluctuation-induced interaction is given by 



V^{R,0i,92) = - 



128 



t2t2 
i?4 



X COS' 



[2 {01 + 02) 



(18) 




FIG. 1. Two rod-shaped inclusions embedded in a membrane. The rods are separated by a distance 7?. 
length Li, width e^, and makes an angle Oi with the line joining the centers of the two rods. 



The ith rod has 



where 0i and 02 are the angles between the rods and the 
line adjoining their centers, as indicated in Fig. |l|. 

The orientational dependence is the square of a 
quadrupole-quadrupole interaction, with the unusual 
property of being minimized for both parallel and perpen- 
dicular orientations of the rods. The above fluctuation- 
induced interactions decay less rapidly at large distances 
than van der Waals forces and may play an important 
role in aligning asymmetric inclusions in biomembranes. 
Since orientational correlations are often easier to mea- 
sure than forces, this result may also be useful as a probe 
of the fluctuation-induced interaction. Finally, this inter- 
action could give rise to novel two-dimensional structures 



for collections of rodlike molecules. In particular, the re- 
semblance of the orientational part of the interaction to 
dipolar forces suggests that a suitable way to minimize 
the energy of a collection of rods is to form them into 
chains. (If the rods are not colinear, the interactions 
cannot be simultaneously minimized.) 

An important property of fluctuation-induced interac- 
tions is that they are non-additive, and cannot be ob- 
tained by adding two-body potentials. For example, con- 
sider an interaction U{\ri ~ r2\)duidu2, between any two 
infinitesimal segments of two rods in Fig. |^. If both rods 
are of length L at a distance L, expanding |ri —r2\, 
and integrating over the two rods, leads to the interaction 
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V{R,9u92) ^ L^U{R) 



(U'{R) 



6 \ R 



+ U"{R)^ (19) 



L* (U'{R) 



12 \ R 



- U"{R) (cos26'i +COS26I2). 



The angular dependence is now completely different, and 
minimized when the two rods are parallel to their axis of 
separation. Presumably both interactions are present for 
rods of finite thickness; the additive interaction is pro- 
portional to L^{L£/ R)^ , where e is the thickness. The 
previously calculated interactions are thus larger by a 
factor proportional to (i?/e)^ and should dominate at 
large separations. 

The case of stiff linear inclusions at close separations 
(L ^ R) is considered in Ref. [Q. It is shown that a 
finite rigidity of the linear inclusions leads to a screen- 
ing out of the Casimir-type fluctuation-induced attrac- 
tion. The screening length is set by the ratio between the 
rigidity of the polymer and that of the membrane. This 
is the length scale below which the polymers are seen 
as straight parallel lines, hence resulting in a Casimir 
interaction. Moreover, the attractive interactions could 
lead to an instability in the shape of the stiff polymers 
signalling a major reduction in their rigidity (softening) 
induced by the membrane fluctuations M] . 



IV. ROUGH SURFACES 

Most computations of Casimir forces are for simple ge- 
ometries, e.g. between two parallel plates, or perfect 
spheres. It is natural to consider how these forces are 
modified by the roughness that is present in most "ran- 
dom" surfaces. There are a number of calculations that 
go beyond the simple planar geometry. For example, in 
Ref. a multiple scattering approach is used to com- 
pute the interactions for arbitrary geometry in a pertur- 
bation series in the curvature. A generalization of the 
approach due to Dzyaloshinskii, Lifshitz, and Pitaevskii 
1^ is developed in Ref. |Q to study the Casimir forces 
for surfaces with roughness. A phenomenological ap- 
proach is introduced in Refs. in which small devi- 
ations from plane parallel geometry are treated by using 
an additive summation of Casimir potentials. However 
(as demonstrated in the previous section) , fluctuation in- 
duced forces are not additive, and additional steps are 
necessary to correct the result p9|]. Another perturbative 
approach is introduced in Ref. ]50||, which could in prin- 
ciple be used to treat surfaces with roughness, although 
not explicitly carried out in this paper. Most of these 
approaches suffer from rather cumbersome treatments of 
the boundary conditions. 

In Ref. Il^, a path integral approach is introduced 
that makes possible relatively simple computations of the 
fluctuation induced force. (A sketch of this formulation 
is presented in Appendix ^) This approach has a num- 
ber of advantages. First, different manifolds (with arbi- 



trary intrinsic and embedding dimensions) in various cor- 
related fluids can be treated in a similar fashion. Second, 
the boundary conditions are quite easily implemented, 
and the corrections can be computed perturbatively in 
the deformations. While this method was originally de- 
veloped for the study of thermal fluctuations, it can be 
adapted to quantum fluctuations, as discussed in the next 
section. In the remainder of this section we calculate the 
corrections to the thermal Casimir force due to substrate 
roughness. 

Many solid surfaces produced by rapid growth or depo- 
sition processes are characterized by self-similar fluctua- 
tions [Q. The fluctuations of a self-affine surface grow 



[h{^)-h{yW^As\^-y 



2Cs 



(20) 



where the overbar denotes quenched average, and is 
a characteristic roughness exponent. The Casimir force 
between a flat and a such a rough surface (with a cor- 
related fluid in between) is calculated in Ref. |15|. The 
resulting free energy per unit area is 



' m IGtt iJ'i 167r ij4-2Cs 4 

(21) 



where Ci is a numerical coefficient [|l5| , and L is the ex- 
tent (upper cutoff) of the self-afhne structure, satisfying 

AH = A'lL'^s < H. (This is the condition that the to- 
tal width due to roughness, AH is less than the average 
separation H , so that the plates are not in contact.) As 
long as L ^ H ^ AH, the interactions in Eq.(^) are 
arranged in order of increasing strength. The largest ef- 
fect of randomness is to increase the Casimir attraction 
by an amount proportional to (AH/H)^. There is also 
another correction term, of the opposite sign, that decays 
as l/_ff'*^^^s, and in principle can be used to indirectly 
measure the roughness exponent ^5. In Eq. (pO|) , if all 
lengths are measured in units of an atomic scale oq (e.g. 
the diameter of a surface atom). As becomes dimension- 
less. Using a reasonable set of parameters: « 0.35, 
ao ~ SA, As ~ 1 and L « 300A, we estimate that for sur- 
faces of 1mm size, and lOOA apart, the forces generated 
by the three terms in Eq.(|lJ) are 1.9 x 10"", 4.9 x 10"^, 
and 3.7 x 10~^N respectively, (using a reasonable lower 
cutoff of ~ 20A), which is measurable with the current 
force apparatus |9|. 

V. THE DYNAMIC CASIMIR EFFECT 

A. Background 

Although less well known than its static counterpart, 
the dynamical Casimir effect, describing the force and 
radiation from moving mirrors has also garnered much 
attention p2-M. This is partly due to connections 
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to Hawking and Unruh effects (radiation from black 
holes and accelerating bodies, respectively) , suggesting a 
deeper link between quantum mechanics, relativity, and 
cosmology [^HJsj. 

The creation of photons by moving mirrors was first ob- 
tained by Moore [£2| for a 1 dimensional cavity. Fulling 
and Davis p3[ demonstrated that there is a correspond- 
ing force even for a single mirror, which depends on the 
third time derivative of its displacement. These computa- 
tions take advantage of conformal symmetries of the 1-1-1 
dimensional space time, and can not be easily generalized 
to higher dimensions. Furthermore, the calculated force 
has causality problems reminiscent of the radiation reac- 
tion forces in classical electron theory [|4j. It has been 
shown that this problem is an artifact of the unphysical 
assumption of perfect reflectivity of the mirror, and can 
be resolved by considering realistic frequency dependent 
reflection and transmission from the mirrors p^ . 

Another approach to the problem starts with the fluc- 
tuations in the force on a single plate. The fluctuation- 
dissipation theorem is then used to obtain the mechanical 
response function [^5| , whose imaginary part is related to 
the dissipation. This method does not have any causality 
problems, and can also be extended to higher dimensions. 
(The force in H-3 dimensional space-time depends on the 
fifth power of the motional frequency.) The emission of 
photons by a perfect cavity, and the observability of this 
energy, has been studied by different approaches [p6|~p8[ . 
The most promising candidate is the resonant production 
of photons when the mirrors vibrate at the optical reso- 
nance frequency of the cavity ||59| ]. A review, and more 
extensive references are found in Ref . | |60| . More recently, 
the radiation due to vacuum fluctuations of a collapsing 
bubble has been proposed pl]-|6^ as a possible explana- 
tion for the intriguing phenomenon of sonoluminescense. 
(Subsequent experimental measurements of the duration 
of the signal may favor more classical explanations.) 

A number of authors have further discussed to notion 
of frictional forces: Using conformal methods in 1+1 di- 
mensions, Ref. |65|l finds a friction term 



(22) 



i^friction(if) = a Fst.ticiH) H/c 



for slowly moving boundaries, where a is a numerical 
constant that only depends on dimensionality. The ad- 
ditional factor of [v/cY' makes detection of such a force 
yet more delicate. There are a few attempts to calculate 
forces (in higher dimensions) for walls that move laterally, 
i.e. parallel to each other [|66|-|70|: It is found that bound- 
aries that are not ideal conductors experience a friction 
as if the plates are moving in a viscous fluid. This friction 
has a complicated dependence on the frequency depen- 
dent resistivity of the plates, and vanishes in the cases 
of ideal (nondissipating) conductors or dielectrics. The 
"dissipation" mechanism for this "friction" is by induc- 
ing eddy currents in the nonideal conductors, and thus 
distinct from the Casimir effect. Possible experimental 



evidence of such a contribution to friction has been re- 
cently reported in Ref. [f7l| . The experiment employs 
a quartz crystal microbalance technique to measure the 
friction associated with sliding of solid nitrogen along a 
lead surface, above and below the superconducting tran- 
sition temperature of lead. An abrupt drop in friction is 
reported at the transition point as the substrate enters 
the superconducting state [ 7l| . 

An interesting analog of the dynamic Casimir effect is 
suggested for the moving interface between two different 
vacuum states of superfluid 3He In this system, the 
Andreev reflection of the masslcss "relativistic" fermions 
which live on the A-phasc of the interface provides the 
corresponding mechanism for friction: The interface is 
analogous to a perfectly reflecting wall moving in the 
quantum vacuum. 



B. Path integral formulation 

The path integral methods originally developed for 
rough surfaces j73j can also be applied to the problem 
of perfectly reflecting mirrors that undergo arbitrary dy- 
namic deformations. Consider the path integral quanti- 
zation of a scalar field 6 with the action 



S 



1 



(23) 



where summation over fi = 0, • • ■ ,3 is implicit. Follow- 
ing a Wick rotation, imaginary time appears as another 
coordinate X/^ ~ ict, in the 4-dimensional space-time. In 
principle, we should use the electromagnetic vector po- 
tential Afj_{X), but requirements of gauge fixing compli- 
cate the calculations, while the final results only change 
by a numerical prefactor. (We have explicitly reproduced 
the known result for gauge fields between flat plates by 
this method Q.) We would like to quantize the field 
subject to the constraints of its vanishing on a set of 
n manifolds (objects) defined by A" = Xa{ya), where 
ya parametrize the ath manifold. We implement the 
constraints using delta functions, and write the partition 
function as 



^ n 

Z = J Vcj,{X)l[l[S{cl>{XM)) exp 



a=l y. 



SI 



(24) 



The delta functions are next represented by integrals 
over Lagrange multiplier fields. Performing the Gaussian 
integrations over (j>{X) then leads to an effective action 
for the Lagrange multipliers which is again Gaussian . 
Evaluating Z is thus reduced to calculating the logarithm 
of the determinant of a kernel. Since the Lagrange mul- 
tipliers are defined on a set of manifolds with nontrivial 
geometry, this calculation is generally complicated. To 
be specific, we focus on two parallel 2d plates embedded 
in 3+1 space-time, and separated by an average distance 
H along the a;3-direction. Deformations of the plates 
are parametrized by the height functions hi (x, t) and 
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/i2(x, t), where x = {xi,X2) denotes the two lateral space 
coordinates, while t is the time variable. As sketched in 
Appendix In Z can be calculated by a perturbative se- 
ries in powers of the height functions . The resulting 



expression for the effective action (after rotating back to 
real time), defined by SoS = —ihlnZ, and eliminating h 
independent terms, is 



-A^{q,uj) (/ii(q, w)/i2(-q, -w) + hi{-q, -uj)h2{(l,uj))] + 0{h^). 



C. The response function 

The kernels A±{q,LLj), that are closely related to the 
mechanical response of the system (see below) , are func- 
tions of the separation H, but depend on q and u only 



through the combination Q ~ q 



7c^ The closed 



forms for these kernels involve cumbersome integrals, and 
are not very illuminating. Instead of exhibiting these 
formulas, we shall describe their behavior in various re- 
gions of the parameter space. In the limit H ^ oo, 
A^(g,w) = 0, and 



ATiq,L0) 



1 



3607r2c5 



-(c^5^ - up-y'l'^ for uj < cq, 

i sgn(a;)(ci;^ — c^q^)^!"^ for cj > cq, 

(26) 



where sgn(w) is the sign function. While the effective 
action is real for > 0, it becomes purely imaginary 
for < 0. The latter signifies dissipation of energy 
p5[ , presumably by generation of photons Q . It agrees 
precisely with the results obtained previously |Q for the 
special case of flat mirrors (q = 0). (Note that dissipa- 
tion is already present for a single mirror.) 

In the presence of a second plate (i.e. for finite H), the 
parameter space of the kernels subdivides into three dif- 
ferent regions as depicted in Fig. 2. In region I {Q^ > 
for any H) , the kernels are finite and real, and hence there 
is no dissipation. In region Ila where — 7r^/iJ^ < < 0, 
the iJ-independent part of is imaginary, while the 
i?-dependent parts of both kernels are real and finite. 
(This is also the case at the boundary = —tt^/H'^.) 
The dissipation in this regime is simply the sum of what 
would have been observed if the individual plates were 
decoupled, and unrelated to the separation H. By con- 
trast, in region lib where < —tt'^/H'^. both kernels 



diverge with infinite real and imaginary parts |75|. This 



i/-dependent divergence extends all the way to the nega- 
tive axis, where it is switched off by a prefactor. 



' Q2=q2- 


co2 

I 






Ila 


>■ 




lib 





FIG. 2. Different regions of the (q, <jj) plane. 

As a concrete example, let us examine the lateral 
vibration of plates with fixed roughness, such as two 
corrugated mirrors. The motion of the plates enters 
through the time dependencies /ii(x, i) = /ii(x — r{t)) 
and ft,2(x, i) — /i2(x); i.e. the first plate undergoes lat- 
eral motion described by r{t), while the second plate is 
stationary. The lateral force exerted on the first plate 
is obtained from fi{t) — SSes/6ri{t). Within linear re- 
sponse, it is given by 



where the "mechanical response tensor" is 



Xi]{uj) = he 



^ q,q, |[A+(q,c.) - A+{q,0)] \h,{q)\^ + ^A^{q,0) (/^i (q)/i2 (-q) + (-q)/i2 (q)) } , 



and there is a residual force 



he 



2TrS{Lu) 



d^q 
(2^ 



%A_(g,0) (/ii(q)/i2(-q) - hi{-q)h2{q)) . 



(27) 



(28) 



(29) 
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VI. CORRUGATED MIRRORS 

Let us now consider a corrugated plate with a de- 
formation /ii(x) — (icos(k-x). From the frequency- 
wavevector dependence of the mechanical response func- 
tion in Eq.(Eq), we extract a plethora of interesting re- 
sults, some of which we discuss next. 



A. Mass corrections 

For a single plate {H — > oo), we can easily calcu- 
late the response tensor using the explicit formulas in 
Eq. (^6|) . In the limit of w <^ cfc, expanding the re- 
sult in powers of to gives Xij = SniijUj'^ + O{oj'^), where 
Sniij — Ah(Pk^kikj/{288Tr^c), can be regarded as correc- 
tions to the mass of the plate. (Cut-off dependent mass 
corrections also appear, as in Ref. |60|.) Note that these 
mass corrections are anisotropic with 



Sm 



1 



2887r2 



-Ad^k 



2;,5 



6m± = 0. 



(30) 



Parallel and perpendicular components are defined with 
respect to k, and A denotes the area of the plates. The 
mass correction is inherently very small: For a macro- 
scopic sample with d k, X = 2'K/k k, 1mm, density 
K, 15gr/cm , and thickness t k, 1mm, we find 5m/ m ^ 
10"'^'*. Even for deformations of a microscopic sample of 
atomic dimensions (close to the limits of the applicabil- 
ity of our continuum representations of the boundaries), 
(5m/m can only be reduced to around 10^^*^. 

With the second plate at a separation H, the mass 
renormalization becomes a function of both k and H , 
with a crossover from the single plate behavior for 
kH ^ 1. In the limit of kH 4C 1, we obtain 5ra\\ — 
hABk^d'^/48cH^ and (5m j_ = 0, with B = -0.453. Com- 
pared to the single plate, there is an enhancement by 
a factor of {kH)~^ in 5m\\. While the actual changes 
in mass are immeasurably small, the hope is that its 
anisotropy may be more accessible, say by comparing 
oscillation frequencies of a plate in two orthogonal di- 
rections. 



B. Dissipation: 

For LD ^ ck the response function is imaginary, and we 
define a frequency dependent effective shear viscosity by 
Xij(w) = —iujriij{u}). This viscosity is also anisotropic, 
with 



1 



7207r2 c4 



-Ad^k 



2 1^2, ,4 



(31) 



Note that the dissipation is proportional to the fifth time 
derivative of displacement, and there is no dissipation for 
a uniformly accelerating plate. However, a freely oscillat- 
ing plate will undergo a damping of its motion. The char- 
acteristic decay time for a plate of mass M is r « 2M/t]. 
For the macroscopic plate of the previous paragraph, vi- 
brating at a frequency of oj « 2ck (in the lO^^Hz range), 
the decay time is enormous, r ~ 10^*s. However, since 
the decay time scales as the fifth power of the dimension, 
it can be reduced to 10~^2s, for plates of order of 10 
atoms. However, the required frequencies in this case (in 
the lO^^Hz range) are very large. Also note that for the 
linearized forms to remain valid in this high frequency 
regime, we must require very small amplitudes, so that 
the typical velocities involved v ~ tquj, are smaller than 
the speed of light. The effective dissipation in region Ila 
of Fig. 2 is simply the sum of those due to individual 
plates, and contains no H dependence. 



C. Resonant Emission: 



The cavity formed between the two plates supports 
a continuous spectrum of normal modes for frequencies 
uj'^ > c2(fc2 -t- TT^/iJ^). We find that both real and imag- 
inary parts of A±{k,Lu), diverge in this regime, which 
we interpret as resonant dissipation due to excitation of 
photons in the cavity. Resonant dissipation has profound 
consequences for motion of plates. It implies that due 
to quantum fluctuations of vacuum, components of mo- 
tion with frequencies in the range of divergences cannot 
be generated by any finite external forcel The imaginary 
parts of the kernels are proportional to the total num- 
ber of excited photons [Q. Exciting these degrees of 
motion must be accompanied by the generation of an in- 
finite number of photons, requiring an infinite amount of 
energy, and thus impossible. However, as pointed out in 
Ref. , the divergence is rounded off by assuming finite 
refiectivity and transmissivity for the mirrors. Hence, in 
practice, the restriction is softened and controlled by the 
degree of ideality of the mirrors in the frequency region 
of interest. 

Related effects have been reported in the literature for 
1+1 dimensions [5^-|59|], but occuring at a discrete set of 
frequencies ojn = mrc/ H ^ with integer n>2. These res- 
onances occur when the frequency of the external pertur- 
bation matches the natural normal modes of the cavity, 
thus exciting quanta of such modes. In one space di- 
mension, such modes are characterized by a discrete set 
of wavevectors that are integer multiples of n/H. The 
restriction to n > 2 is a consequence of quantum electro- 
dynamics being a 'free' theory (quadratic action): only 
two-photon states can be excited subject to conservation 
of energy. Thus the sum of the frequencies of the two 
photons should add up to the external frequency 
In higher dimensions, the appropriate parameter is the 
combination w^/c^ — . From the perspective of the ex- 
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cited photons, conservation of momentum requires that 
their two momenta add up to q, while energy conser- 
vation restricts the sum of their frequencies to lu. The 
in-plane momentum q, introduces a continuous degree of 
freedom: the resonance condition can now be satisfied 
for a continuous spectrum, in analogy with optical res- 
onators. In Ref. [Q, the lowest resonance frequency is 
found to be 2nc/H which seems to contradict our pre- 
diction. However, the absence of oji = ttc/H in 1+1 D 
is due to a vanishing prefactor , which is also present 
in our calculations. However, in exploring the continu- 
ous frequency spectrum in higher dimensions, this single 
point is easily bypassed, and there is a divergence for 
all frequencies satisfying up' jt? > (7^ + 7r^/i?^, where the 
inequality holds in its strict sense. 



metric function with respect to wq/Z, where it is peaked. 
The peak sharpens as the parameter wo/Zcq — > 0, and 
saturates for fco = 0. 

The connection between the dissipative dynamic 
Casimir force, and radiation of photons, is made explicit 
by calculating the total number of photons radiated per 
unit time and per unit area of the plate. The result is 
identical to the energy dissipation rate calculated in Ref. 
[ |73[ . No radiation is observed at frequencies higher than 
tiio I due to conservation of energy, and also for ^0/^0 < 1 : 
in agreement with section B above, where no dissipative 
forces are found in this regime. 
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D. Radiation spectra 

Where does the energy go when the plates experience 
viscous dissipation? When the viscosity is a result of 
losses in the dielectric boundaries ||6^-|70||, the energy is 
used up in heating the plates. Since we have examined 
perfect mirrors, the dissipated energy can only be ac- 
counted for by the emission of photons into the cavity. 
The path integral methods can be further exploited to 
calculate the spectrum of the emitted radiation ||7^ . The 
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tons (at t -l-oo), to two point correlation functions of 
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deformations. From the transition amplitude (after in- 
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Specifically, calculations of the angular distribution 
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perfectly reflecting plate, undulating harmonically with 
a frequency loq^ and a wavevector kg. Depending on the 
ratio wo/fco: we find that radiation at a frequency is 
restricted to a particular window in solid angle . The 
total spectrum of radiation is found by integrating the 
angular distribution over the unit sphere, and is a sym- 
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APPENDIX A: PATH-INTEGRAL 
FORMULATION AND DEFORMED SURFACES 

In this Appendix, we sketch the path integral formu- 
lation developed in Ref. ||l5| for surfaces with roughness. 
Consider n manifolds embedded in a d-dimensional cor- 
related fiuid, with an energy cost appropriately general- 
ized from Eq.(^). The manifolds are described by the 
functions Ra{xa), where Z)Q.-dimensional inter- 

nal coordinate (Da — 1 for a polymer, and Da — 2 for 
a membrane), while Ra indicates a position in the d- 
dimensional fluid. The fluctuation-induced interactions 
between the manifolds are obtained by integrating over 
all configurations of the field 0, subject to the constraints 
of its vanishing on the external manifolds. The bound- 
ary conditions {<i>(Ra{xa)) = 0, for a = l,2,...,n) are 
imposed by inserting delta functions. Using the integral 
representation of the delta function, we can write 



exp 



U 



knT 



= -Y /^<^('') n ^V'a(a;a)exp ^-7io[0] +i /"dx^ 

^ a— 1 ^ 



V'a(a;a)(/)(i?Q (Xa)) 



(1) 



where ipaixa) are the auxiliary fields defined on the n manifolds, acting as sources coupled to </>. For the quadratic 
Hamiltonian of Eq.(^), it is easy to integrate over the field (p, and obtain the long-range interactions between the 
sources as 



exp 



/n 
VlpaiXa) exp {-Hl[tpa{Xa)]) 



The action Ti-ilipaixa)] is a quadratic form for the n component field ^ = {ipi, ip2T ■ ■ ■ given by 

n n „ 



(2) 



(3) 



a = l 13=1 ' 
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where G'^{r) = (<^(r)(/)(0))Q is the two-point correlation 
function of the field (j) in free space. Finally, the effec- 
tive interaction between the manifolds is obtained by per- 
forming the Gaussian integrations over the field ^ as 



Heff[RaiXa) 



\iidet{M[Ra{xa)]) 



(4) 



The matrix M (which can be read off from Eq.(|^)) 
is a functional of Ra{xa) and its determinant is in gen- 
eral difficult to evaluate for arbitrary configurations, ft 
is possible, however, to perturbatively calculate the cor- 
rections due to small deformations around simple base 
configurations. Consider two surfaces in ci = 3, with 
average separation H. One plate has small deforma- 
tions described by /i(x), i.e. -Ri(x) = (xi, 2:2,0), while 



i?2(x) = {xi,X2tH + h{x)). The effective Hamiltonian 
can now be written as Ti-eff = 'Hfiat+'Hcom where Hfiat 
is the Casimir interaction for two flat plates, and Ticorr 
is the additional cost of deformations. For flat plates of 
area A, the interaction energy is 



Ti-flat 

A 



= kBT 



(27r)2 '^\2ap 



C(3) keT 
167r i/2 



(5) 



The first term in Eq.(|^) is a contribution to the surface 
tension which depends on a lattice cutoff. The second 
term is the usual Casimir interaction, decaying as 
and with a universal amplitude — C(3)/167r w —0.02391. 
The energy cost of deformations is given by 



He 



= -knT X 



16 



^ J d^x/^2(x) + d^xdV[Mx) - My)]^ X 



B7r2|x-y|6 



27r|x-y|3ij3 



H 



1 



H 



1 



|x-y| 

H 



(6) 



with two kernel functions defined by 

K,{t) ^!^^u^{^-~l)-^Mtu), (7) 

K2{t) = f u2gu(g2„ _ 1)-1 Jo(to). (8) 

The first term in Eq.(|^) represents an instability to de- 
formations that is related to the attraction between the 
plates. Remarkably, this term can be obtained intuitively 
by replacing l/iJ^ .^^ith l/(i? + /i(x))2 in Eq.(||) and 
averaging over the position x. The second term repre- 
sents long-range interactions between the deformations, 
induced by the fluctuations of the field. The first term in 
the curly brackets is the conformation energy of the de- 
formed surface in the absence of the second plane, and is 
independent of H. The remaining terms represent corre- 
lations due to the presence of second plane. Both Ki{t) 
and K2{t) approach a constant as t ^ 0. As t ^ 00, 
Ki{t) ~ l/<3, while K2{t) - exp(-M), with h « 3.3. 
The large t behaviors of Ki (t) and K2 (t) determine the 
long-range interactions between height fluctuations. 

APPENDIX II: PATH INTEGRAL 
FORMULATION OF CHARGED FLUIDS 

Here, we introduce a systematic path integral for- 
mulation to study fluctuation-induced interactions in a 
charged fluid. Consider n charged manifolds embedded 
in a d-dimensional aqueous solution of neutralizing coun- 
terions, interacting through Coulomb potentials. The 
manifolds have charge densities — Uq (all assumed to be 



negatively charged for simplicity), and are described by 
the functions Ra{xa), where Xa is a -dimensional in- 
ternal coordinate, while Ra indicates a position in the 
d-dimensional solution. There are Nc positively charged 
counterions of valence z, each described by a position 
vector Ri, in the d-dimensional solution. The Coulomb 
Hamiltonian can be written as 

Hc^lJ d'^Xd'^X' p(X)-^^^^p(X'), (1) 

where 

p{X) = - ^ / dxo.cJc.5\X - Rc^ixa.)) (2) 

+ Y,zS^iX-R^), 
i=l 

is the number density of the charges. Charge neutral- 
ity requires — '^aAa + zNc = 0, where Aa is the 
Da-dimensional area of the ath manifold. 

A restricted partition function of the Coulomb system, 
depending upon the shapes and locations of the macori- 
ons, is now given by 

2nARM]= [U^^""^^''^^ (3) 

i—l 

in which a is a short-distance cut-off. Using the Hubbard- 
Stratanovich transformation of the Coulomb interaction, 

(4) 



e-"c/ksT^ J p^(^)cxp|-|g J d'^X(V0)^+z J d''Xp{X)cf>{X)Y 
we can rewrite the partition function as 
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•^'^ e^^'f'^A''\ (5) 



where Sd is the area of the d-dimensional unit sphere. We can introduce a fugacity y, and a rescaled partition function 
that can be rewritten as 

^ = X/ ^N.N,j;^2:N[Ra{Xa)] 



A shift in the field (f) by ~9, and use of the neutrahty condition renders the ^-integration trivial. We can then sum 
up the exponential series, and obtain 

Z[RM] = J V(P{X) e-^[^l, (8) 



in which 



n[<t>] = ^ J d'*X(V^)2 +iJ2jdx^ <ra<i^{RM) -^J d-^^e-^f^). (9) 

CK — 1 

Note that the fugacity y can be eliminated using the identity 

91nZ , , 

iVe^y^, (10) 

which follows from Eq.(|^). 

We next evaluate the path integral using a saddle point approximation. The extremum of Eq. (^) , obtained from 
5T-L/54) = 0, is the solution of the Poisson-Boltzmann (PB) equation 

- V2(zV(^)) - ^^^e-^'^^^) --ill Ax^^4^^'^^ - ^"(^"))' (11) 

for the (real) field 'ip{X) — —i(f>{X), in which — S^e^yz^ /eksTa'^ defines the inverse square of the Debye screening 
length. To study the fluctuations on top of this saddle point, we can set 4> = (f> + S(j), and expand the Hamiltonian up 
to quadratic order, to get 

= + / d^X [(V<50)2 + «2e-^'^(^)<5(/.2l . (12) 

25^6^ J I -I 



The free energy of the system of charged manifolds in 
the presence of fluctuating counterions now reads 



Ft 



PB 



In det 



(13) 



where FpB = H[z'0(X)] is the Poisson-Boltzmann free 
energy, and m?{X) = K'^e~^^^'^'> is a "mass (or charge 
compressibility) profile" . The PB free energy is known to 
be generically repulsive 26|. The fluctuation- induced 
correction, however, is attractive. For highly charged 
manifolds, it is indeed reminiscent of the Casimir interac- 



tions, but with the boundary constraints smoothed out. 
To see this, one should note that the mass profile is indeed 
identical to the density profile of the counterions. Highly 
charged manifolds accumulate counterions in their vicin- 
ity, and consequently the fluctuations of the "potential" 
field (/> are suppressed in a region close to the manifolds, 
but are unconstrained in other regions in the solution; 
hence leading to a Casimir-type fluctuation-induced at- 
traction. 
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